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ABSTRACT 

We present spectral and long-timescale variability analyses of Chandra X-ray Observatory ACIS- 
S observations of the 19 millisecond pulsars (MSPs) with precisely known positions in the globular 
cluster 47 Tucanae. The X-ray emission of the majority of these MSPs is well described by a thermal 
(blackbody or neutron star hydrogen atmosphere) spectrum with a temperature Toff ~ (1 — 3) x 10^ 
K, emission radius i?off ~ 0.1 — 3 km, and luminosity Lx ^ I0'^''~3i gj-gg g-i_ p^j. several MSPs, there 
is indication that a second thermal component is required, similar to what is seen in some nearby 
field MSPs. The observed radiation most likely originates from the heated magnetic polar caps of the 
MSPs. The small apparent scatter in Lx is consistent with thermal emission from the polar caps of 
a global dipole field although the small emission areas may imply either a more complex small-scale 
magnetic field configuration near the neutron star surface or non-uniform polar cap heating. The radio- 
eclipsing binary MSPs 47 Tuc J, O, and W show a significant non-thermal (power-law) component, 
with spectral photon index F '--^ 1 — 1.5, which most likely originates in an intrabinary shock formed 
due to interaction between the relativistic pulsar wind and matter from the stellar companion. We 
re-examine the X-ray-spindown luminosity relation {Lx — E) and find that for the MSPs with thermal 
spectra Lx oc E^ , where /3 ~ 0.2 ± 1.1. Due to the large uncertainties in both parameters the result 
is consistent with both the linear Lx — E relation and the flatter Lx oc E^'^ predicted by polar cap 
heating models. In terms of X-ray properties, we find no clear systematic differences between MSPs 
in globular clusters and in the field of the Galaxy. We discuss the implications of these results of the 
present understanding of the X-ray emission properties of MSPs. 

Subject headings: globular clusters: general — globular clusters: individual (47 Tucanae) — pulsars: 
general — stars: neutron — X-rays: stars 



1. INTRODUCTION 

Millisecond pulsars (MSPs) represent a distinct popu- 
lation of rotation-powered pulsars, characterized by short 
spin periods, P < 25 ms, and small intrinsic spin-down 
rates. Pi ^ 10"^*^, implying relatively low surface mag- 
netic dipole field strengths Bsurf oc {PPi)^'"^ - 10^"^° G 
and large characteristic spindown ages r = P/2Pi > 1 
Gyr. These objects have been studied e xtensively at 
radio wavelengths since their discovery IjBacker et al.l 
Il982|l but at X-ray energies, where the bulk of ob- 
servable radiation is expected, this has become possi- 
ble only relatively recently, owing to their intrinsic faint- 
ness {Lx < 10^^ ergs s"^). The ROSAT observatory 
was the first to detect X-ray emission from MSPs, in- 
cluding 10 in the field of the Galaxy and one in the 
globular cluster M28 (see Becker & Triimper 1999, and 
references therein). In recent years, the unprecedented 
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sensitivity and spatial resolution of the Chandra X- 
ray Observatory have allowed the detection of X-ray 
counterparts to radio MSPs in several globular clus - 
ters: 47 Tucanae a nd NGC 6397 llGrindlav et alJl2002tl 
NGC 6752 (D'A mico et all l2002[)_, M 4 llBassaetal] 
I2004|) . possibly M30 fRansom_et al 2004), and M71 (Ei- 
sner et al., in preparation). In addition, the Chan- 
dra and XMM-Newton telescopes have allowed much 
more detailed studies of the existing sam ple of X-ray- 
detect ed MSPs than previously possible llZavlin et alJ 
l200^; Becker & Aschenbach 2002: 'Becker et al.' '2003t 
Webb. Olive, & Barret 2004a: Webb et al. 2004b) . 

The X-ray emission detected from MSPs is observed 
to be both of non-thermal and thermal character. In the 
X-ray bright MSPs, B1821-24 and B1937+21, the nar- 
row X-ray pulse profiles indicate the emission is beamed 
and, therefore, must be non-thermal radiation origi- 
natin g in the pulsar magnetospherc (Becker & Tri impeJ 
I1999D . Alternatively, as seen in the "black widow" pul- 
sar, B1957-I-20, non-pulsed non-thermal X-rays can be 
produced by interaction of the relativistic particle wind 
from the pulsar with matter fr om a stellar compan- 
ion or the interstellar medium IjStappers et al. 2003). 
Thermal em i ssion observed from nearby MSPs (sec e.g. 
iZavlin et all l20(il iBecker fc AschenbachI l200l . is be- 
lieved to originate from the magnetic polar caps of the 
underlying neutron star (NS), heated by a backflow of en- 
ergetic particles from the pulsar magnetosphere (se e e.g. 
IHarding fc Muslimovl 12002 IZhang fc Ch^ 12001 and 
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Fig. 1. — A merged 0.3 — 6 keV image of all four Chandra ACIS-S observations of 47 Tucanae showing the positions of the X-ray 
counterparts of 19 millisecond pulsars. The 1" circles are centered on the radio MSP locations and enclose > 90% of the total energy for 
0.5 — 1 keV thermal sources. The grayscale corresponds to number of counts increasing logarithmically from zero (white) to 8306 (black). 
North is up and cast is to the left. 



references therein). 

The detection of 16 MSPs in the globular cluster 47 
Tucanae (NGC 104, henceforth 47 Tuc) by the Chan- 
dra ACIS-I instrument has granted an excellent op- 
portunity for a detailed study of the physical proper- 
ties of this class of objects (Grindlav et al. 2001. 2002) . 
These pulsars we r e discovered a,t radio wavelengths b y 
iManch es ter et al 



and Caniilo et al 



3 aiscoverea at raaio waveien gtns b y 

fT99n[ dp), IRobinsonet'aD IT99I . 

. With furth er pulse timing ob- 



servations. lFreire et al.l f 2001a". "2003"^ were able to obtain 
precise positions for 16 of the 22 MSPs presently known 
in th is globular clust er (Caniilo & Rasio 2005). In addi- 
tion, lEHrnQn^^OU <|2002i) detected the optical counter- 
part of 47 Tuc W, which does not have an accurate radio 
timing position, by matching the binary phase and pe- 
riod of an optical candidate with that of the radio pulsar. 
Finally, timing solutions have recently been obtained for 
47 Tuc R and Y (Freire et al., in preparation) thus in- 
creasing the numb er of MSPs with know n positions in 
this cluster to 19. iGrindlav et alJ (|2002) have reported 
in detail on the analysis of the X-ray emission properties 
of the 47 Tuc MSPs, using the Chandra ACIS-I obser- 
vations mentioned above. However, due to the inherent 



faintness of these sources {Lx 



10 



30-31 



ergs s ^) and 



the limited exposure time (72 kiloseconds of total data 



obtained resulting in 1-30 detected counts for each MSP), 
the exact nature of the X-ray emission could not be deter- 
mined with great confidence, although some interesting 
constraints could be placed. In the meantime, additional 
X-ray observations of 47 Tuc have been carried out, al- 
lowing further study of the properties of globular cluster 
MSPs as well as MSPs in general. 

In this paper, we present spectral and long-timcscale 
variability analyses of the X-ray counterparts of MSPs in 
the globular cluster 47 Tuc based on deep Chandra ACIS- 
S obse rvations. This p aper extends the work presented 
in Hci nke et al.l l)2005j) . The work is organized in the 
following manner: §2 describes the observations and data 
reduction procedure; in §3 we discuss the X-ray spectral 
properties of the MSPs, while in §4 we investigate the 
long-term temporal behavior of the X-ray emission; in 
§5 we compare the X-ray and radio pulsar properties. 
Finally, in §6 we discuss the implications of the results 
presented and offer conclusions in §7. 



2. OBSERVATIONS AND DATA REDUCTION 
The data set presented here consists of four 



65- 



kilosccond observations of the core of the globular cluster 
47 Tuc with the S3 chip of the Chandra ACIS-S detector 
at the focus. The first three observations were performed 
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between 2002 September 29 and October 3, while the 
fourth was carried out on 2002 October 11. These ob- 
servations have resulted in the detection of ~300 X-ray 
sources within the half-mass radius of 47 Tuc, including 
X-ray counterparts to all of the 19 MSPs with known 
positions. In this present paper we focus on the counter- 
parts of the radio MSPs; for a comprehensive ana lysis of 
the other X-ray sources, see iHeinke et all ll2005j) . The 
initial data reduction and image processing were per- 
formed using the CIAO^ 3.0 software package and are 
also described in detail bv Heinke et al.i ^2005). Figure 1 
shows the combined image of all four observations with 
1" circles centered at the radio position of each MSP. The 
CIAO tool WAVDETECT was able to detect an X-ray 
source at the position of each of the radio MSPs, includ- 
ing 47 Tuc C, L, M, Q, and T, which were not detected 
or onl y marginally det ected in the 2000 ACIS-I observa- 
tions llCrindlav et"ani200^ . The MSPs 47 Tuc G and 
I, separated by just 0'.'12, could not be individually re- 
solved, whereas 47 Tuc F and S, separated by 0'.'7, were 
again identified as a single elongated source. In this pa- 
per, we also include analy ses of MSP W, whose optical 
counterpart was found bv lEdmonds et all i2002D to co - 
incide with the X-ray source W29 IjGrindlav et al.ll200H) . 
and the MSPs R and Y, whose recently obtained radio 
timing solutions (Freire et al. , in preparation) place them 
at the positions of the X-ray sources W198 and W82, re- 
spectively. 

To facilitate the spectral and variability analyses, the 
counts for each MSP in the 0.3-8 keV energy band were 
extracted from circular regions with 1" radii around the 
radio MSP positions. These circles enclose >90% of the 
total energy for 0.5-1 keV thermal sources. A slightly 
more involved extraction method was required for 47 Tuc 
F and S, whose extraction circles overlap, as evident in 
Figure 1. The relative contribution of counts from each 
MSP was determined by maximum likelihood fits of two- 
dimensional Gaussians to the five X-ray sources in the 
vicinity of the radio positions of F and S using the Sherpa 
application in CIAO. This procedure was carried out in 
the 0.3-4 kcV range, as well as in three narrow bands 
(0.3-0.8 keV, 0.8-1.5 keV, and 1.5-4 keV) to allow an 
approximate spectral analysis. Using this method we 
find that 62 ± 7% of the total counts can be attributed 
to S. For 47 Tuc G and I, this was not possible as their 
radio positions are only 0'.'12 apart. However, the count 
rate observed from their position is roughly twice that 
of a typical MSP in the sample and is, thus, consistent 
with being combined emisson from both MSPs. The net 
counts for most MSPs were obtained by subtracting the 
background rate, obtained from three source-free regions 
on the image, from the total. 

In general, the background noise in the energy range of 
interest was found to be very low (~1 background count 
per pixel over ^^260 ks of total exposure time). This was 
not the case for 47 Tuc L, which is located ver y near 
the bright X-ray source X9 I'Grindl av eran i20011 The 
background for this MSP was obtained by considering a 
section of an annulus of width 2" around the positions 
of X9, at the distance to L, while excluding the 1" ex- 
traction circle around the MSP. Using this background 
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cts(1. 5-4.0 keV)/cts(0.6-1.5 keV) 



^ Chandra Interactive Analysis 
I http://asc.harvard.edu/ciao/t 



of Observations 



Fig. 2. — X-ray color-color diagram for the MSPs in 47 Tuc. 
The lines represent model tracks for absorbed power-law (dotted), 
blackbody (dashed), and neutron star hydrogen atmosphere (dot- 
dashed) spectra convolved with the ACIS-S instrument response. 
The open symbols mark values for the model tracks in intervals of 
0.025 keV for the blackbody (squares) and hydrogen atmosphere 
(circles) temperatures and 0.25 for the power-law photon index 
(triangles). 



region we find that ~30% of the gross number of counts 
are due to X9. For 47 Tuc R, situated 1'.'3 from the 
bright source W24, this procedure was not possible due 
to a combination of the smaller distance to W24 and the 
large pixel size of the ACIS detector. However, the dis- 
tribution of counts within the 1" extraction region of R 
suggests that W24 contributes with no more than 10% 
to the total counts. In the case of MSP O, the degree of 
contamination is difficult to gauge due to severe crowd- 
ing. In addition, the X-ray counterpart of 47 Tuc O, 
W39, is very likely a blend of the MSP and a variable 
source (see §4). 

3. SPECTRAL ANALYSIS 

To get a general sense of the X-ray properties of 
the MSPs, we have computed the hardness ratios 
HR1= (medium counts) /(soft counts) and HR2=(hard 
counts) /(medium counts) for each object and con- 
structed an X-ray color-color diagram. The soft (0.3-0.6 
keV), medium (0.6-1.5 keV), and hard (1.5-4 keV) en- 
ergy bands were selected so as to reduce the errors in 
HR2, allowing us to discriminate between the two fa- 
vored spectral models: thermal or non-thermal (power- 
law). The net counts in the three bands for each MSP 
are listed in Table 1, while the resulting X-ray color-color 
diagram is shown in Figure 2. For comparison, we have 
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TABLE 1 

Net counts from the 19 MSPs in 47 Tuc. 



MSP 



Soft 
(0.3-0.6 keV) 



Medium 
(0.6-1.5 keV) 



Hard 

(1.5-4 keV) 



Total Cnts 
(0.3-4 keV) 



c 


12 ±4 


29 ±5 


4±2 


45 ±7 


D 


20 ±5 


60 ±8 


12 ±4 


92 ±10 


E 


39 ±6 


74 ±9 


11±3 


124 ±11 


F-t-S^" 


67 ±8 


159 ± 13 


20 ±5 


246 ± 16 


pa 


15 ±12 


83 ±21 


12 ±7 


93 ±18 


G-HI 


40 ±6 


102 ± 10 


16 ±4 


158 ± 12 


H 


24 ±5 


46 ±7 


2±2 


72 ±9 


J 


25 ±5 


80 ±9 


32 ±6 


137 ±12 


L 


37 ±8 


132 ± 14 


11 ±7 


180 ± 18 


M 


7±3 


47 ±7 


2±2 


56 ±8 


N 


20 ±5 


40 ±7 


10 ±3 


70 ±9 


O 


30 ±6 


61 ±9 


10 ±4 


101 ± 12 


Q 


18 ±4 


40 ±6 


9±3 


67 ±8 


R 


41 ±8 


105 ±11 


16 ±5 


162 ± 14 


S"" 


52 ±12 


76 ±21 


8±7 


153 ± 20 


T 


10 ±3 


22 ±5 


6±3 


38 ±7 


U 


10 ±3 


67 ±8 


11±3 


88 ±9 


W 


60 ±8 


160 ± 13 


86 ±9 


306 ± 18 


Y 


16 ±5 


41 ±7 


5±3 


62 ±9 



^For F and S, we quote the combined counts as well as the 
counts for each MSP apportioned using the procedure de- 
scribed in the text. 



also plotted values of HRl and HR2 for simple absorbed 
models of blackbody, neutron star hydrogen atmosphere 
(Oioman i 1987 : Zavlin. Pavlov, fc Shibanov 1996; Llovd 
l2003j) 7and power-law spectra for a range of temperatures 



and photon indices, respectively. 

The Chandra ACIS-I observations of 47 Tuc revealed 
that the X-ray colors of the majority of the MSPs are con- 
sistent with those of soft, presumably thermal, sources 
with blackbody temperatures of kT ^ 0.22 keV an d 



emission radii of _R < 0.6 km (jGrindla 
The X-ray color-color diagram for the ACIS-S data 
(Fig. 2) shows that, with the exception of 47 Tuc 
J and W, the MSPs are indeed soft sources although 
the positions of the MSPs on the diagram are sugges- 
tive of a more complex spectrum, such as a compos- 
ite thermal and non-thermal (power-law) spectrum or a 
multi-temperature thermal spectrum. As'Gri ndlav et id] 
1^02) and Ransom et al. (2004) argue, although consis- 
tent with the X-ray colors of the MSPs, a pure thermal 
bremsstrahlung spectrum is ruled out as it would require 
unrealistically large plasma densities, which would be in- 
consistent with the observed dispersion measures. 

The fivefold improvement in counting statistics, result- 
ing from the longer exposure time and the better soft 
response of the ACIS-S compared to the ACIS-I obser- 
vations, permits a more precise determination of the na- 
ture of the X-ray emission via spectral fitting for the 47 
Tuc MSPs. For this purpose, we used the total counts 
extracted for each MSP from the four observations to 
generate source and background spectra, along with the 
corresponding response matrix and auxiliary files. The 
MSP pairs 47 Tuc F and S, and G and I, were treated 
as single sources as the individual spectra could not be 
separated. For the fainter MSPs, the extracted counts 
were grouped into energy bins containing at least 10 
counts, while for the brighter MSPs each bin contained 



at least 15 counts. The subsequent least-squares fits to 
the spectra were performed using the XSPEC^ software 
package and were restricted to photon energies between 
0.3 and 8 keV. In the fitting process, we considered the 
following spectral models as plausible physical descrip- 
tions of the MSP emission: a pure blackbody (BB), an 
unmagnetiz ed neutron star hydrogen atmosphere (NSA; 
lLlovdll200^ . a pure power-law (PL), a composite ther- 
mal and power-law (BB-I-PL or NSA-I-PL), and a two- 
temperature thermal model (BB-I-BB or NSA-I-NSA). 
Throughout this analysis the hydrogen column density 
was fixed at the nominal value for 47 Tuc, iVn = (1.3 ib 
" llCra,tton et a,l " " 



. . - '(-^Tn— IH4ra,tton et alJ-^iMK^predehl fc SchTuittl 

IT995t ICardelli. Clavton. fc Mating I1989D . 



NSA model we 



(2GAfjvs) / [c'Rns 



assumed 

1/2 



a value of z„ 



0.3) X 10^° I 

For the 

'9 - [1 - 
1 — 0.31 for the gravi- 
tational redshift at the NS surface, appropriate for a 
Mns = 1-4 Mq, i?NS = 10 km star. Thus, the only 
free parameters for the BB and NSA fits were the effec- 
tive temperature and BB/NSA normalization, while for 
the PL model the free parameters were the photon index 
and the normalization. 

From the spectral fits we find that many of the 47 
Tuc MSP spectra are well described by a purely thermal 
spectrum (BB or NSA) . For the BB model the fits yielded 
effective temperatures and radii of Tcff ^ (1.5 — 3.0) x 10^ 
K and i?off ^ 0.1 — 0.5 km, respectively, while in the 
case of the NSA model we obtained Toff ~ (0.8 — 1.8) x 
10^ K and i?cff ~ 0.3 - 3 km. Tables 2 and 3 list the 
best fit model, the corresponding values for the fitted 
parameters, and the unabsorbed fiux for each MSP for 
the BB and NSA model fits, respectively. Figure 3 shows 
representative thermal spectra for the 47 Tuc MSPs. 

It is apparent from the values of xt listed in Tables 
2 and 3, as well as Figure 3, that the spectra of sev- 
eral MSPs are poorly fit by a single temperature model, 
especially at energies above ^--^1.5 keV. In addition, we 
find the fits to be unacceptable for a PL or a compos- 
ite PL-t-BB/PL+NSA model as well. Moti vated by the 
result s for the nearby MSPs J0437-4715 (iZaylin et alJ 
1200^ and J0030+0451 IjBecker fc Aschenbachll2002() we 
attempted fitting a two-temperature polar cap model to 
the brighter thermal-emitting 47 Tuc MSPs, namely, D, 
E, L, and R. We found that a two-temperature spec- 
trum can account for the excess at >1.5 keV apparent 
in Figure 3, resulting in an improved fit, although the 
limited count statistics do not allow us to reliably con- 
strain the emission areas. Nonetheless, if we fix the two 
temperatures to values comparable to those of J0437- 
4715 and J00304-0451 we obtain very similar results for 
the effective radii of the 47 Tuc MSPs to those obtained 
for the two field MSPs. Specifically, in the case of the 
NSA model, for Tj= 2.1 x 10^ K and T2 = 0.54 x 10^ 
K (Zavlin et al. 2oS2) the fits yielded Ri ~ 0.3 km and 
i?2 ~ 2 km, respectively, while for a BB model with T-\ = 
3.2 X 10^ K and T2 = 1.4 x 10^ K l)Becker fc AschenbachI 
120021) ■ we obtain Ri - 50 m and R2 ~ 300 m. If we 
extend the analogy further, it seems plausible that the 
thermal emission from all 47 Tuc MSPs consists of two 
components as well. This imples that the derived values 
of Toff and RcB listed in Tables 2 and 3 for the single 

* http://heasarc.gsfc.nasa.gov/docs/xanadu/xspec/ 
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TABLE 2 
Best fit spectral models and unabsorbed fluxes for the 47 Tuc MSPs. 





Spectral 


TrM 


Res'" 


Photon 


Xi/doi 


Fx (0. 


3-8 keV) 


MSP 


Model'' 


(lOf- K) 


(km) 


Index 




(10-15 ergs cm-2 g-i) 


C 


BB 


2.02 + 0.18 


0.11 + 0.08 


- 


0.29/3 


f. 7-1+0.14 
U-' 1-0.26 




D 


BB 


2.20 + 0.17 


0.13 + 0.06 


- 


1.96/8 


i-37t°:J« 




E 


BB 


1.75 + 0.09 


0.28 + 0.17 


- 


1.60/6 


2.08l°:fe 




F+S 


BB 


2.19 + 0.09 


0.22 + 0.11 


- 


1.05/10 


Q on+0-31 




G+I 


BB 


2.21 + 0.12 


0.18 + 0.10 


- 


0.96/12 


9 c /1 0.19+ 




H 


BB 


1.93 + 0.13 


0.17 + 0.11 


- 


0.55/6 


1 so+°i2 




J 


BB+PL 


1.73 + 0.21 


0.22 + 0.17 


1.00 + 0.56 


1.07/5 


4 77+0.95 
*■' '-1.51 


(3.28to-9) 


L" 


BB 


2.27 + 0.10 


0.20 + 0.10 


- 


2.04/11 


o c^+0.32 




M 


BB 


2.22 + 0.18 


0.11 + 0.07 


- 


0.40/4 


1 m+'''i'' 

l-Ul-0.30 




N 


BB 


2.07 + 0.18 


0.13 + 0.09 


- 


1.94/5 


n Qs+0-1'^ 




o= 


BB+PL 


1.76 + 0.15 


0.28 + 0.18 


1.33 + 0.79 


1.21/10 


. /I/I+1.44 
4-44_Qgg 


(2.231?:?^) 


Q 


BB 


2.24 + 0.20 


0.11 + 0.07 


- 


1.17/5 


1 no+°i3 




R 


BB 


2.51 + 0.16 


0.15 + 0.08 


- 


2.54/7 


r, oy + 0.17 




T 


BB 


1.56 + 0.16 


0.19 + 0.17 


- 


2.09/2 


0.63+111 




U 


BB 


2.76 + 0.22 


0.08 + 0.06 


- 


0.55/6 


1 00+0.12 
1-^^-0.33 




W 


BB+PL 


1.52 + 0.28 


0.29 + 0.29 


1.36 + 0.24 


1.22/14 


10.90+0*; 


(9.55lMi) 


Y 


BB 


2.24 + 0.18 


0.11 + 0.07 


- 


1.16/4 


1 n3+''io 

l-'J'^-0.28 





^PL is a power-law and BB is a single temperature blackbody model. For all MSPs the hydrogen column density was fixed 
at A'^H = (1.3 ± 0.3) X 10^" cm^^. In the last column, the values in parentheses for J, O, and W represent the flux in the PL 
component. All uncertainties quoted are la. 



TABLE 3 

Best fit spectral models and unabsorbed fluxes for the 47 Tuc MSPs. 



MSP 



Spectral 

Model'' 



Tctt 
(10*' K) 



RcH 
(km) 



Photon 
Index 



x2/dof 



Fx (0.3-8 keV) 



C 
D 

E 

F+S 

G+I 

H 

J 

L'^ 

M 

N 

o= 
Q 

R 

T 
U 

W 
Y 



NSA 

NSA 

NSA 

NSA 

NSA 

NSA 
NSA+PL 

NSA 

NSA 

NSA 
NSA+PL 

NSA 

NSA 

NSA 

NSA 
NSA+PL 

NSA 



1.12 
1.29 
0.88 
1.27 
1.28 
1.04 
0.89 
1.42 
1.27; 
1.20 
0.98 
1.30 
1.54 
0.80 
1.82 
0.94 
1.35 



0.23 
0.19 
0.20 
0.17 

+0.12 
0.11 

+0.10 
0.09 

+0.12 
0.12 
0.16 
0.14 
0.20 
0.16 

+0.12 
0.12 

+0.19 
0.17 
0.22 
0.19 

+0.17 
0.15 

+0.26 
0.21 

+0.19 
0.17 

+0.19 
0.15 

+0.27 
0.24 

+0.34 
0.25 

+0.20 
0.17 



0.57 + 0.56 
0.61 + 0.51 
1.75 + 1.41 
1.03 + 0.64 
0.84 + 0.58 
0.94 + 0.79 
1.43 + 1.40 
0.78 + 0.50 
0.53 + 0.45 
0.61 + 0.57 
1.48 + 1.27 
0.51 + 0.46 
0.60 + 0.45 

0-2.4 
0.28 + 0.22 

0-2.1 
0.47 + 0.40 



1.0 



1.3 



1.15 + 0.34 



0.68/3 
1.40/8 
1.03/6 
0.77/10 
0.60/12 
0.59/6 
1.00/6 
2.09/11 
0.45/4 
1.38/5 
1.34/11 
0.78/5 
1.60/7 
1.52/2 
0.49/6 
1.17/14 
0.98/4 



0.73 

1.54 

2.28 

4.09 

2.78 

1.39 

4.69 

3.70 

1.08 

1.09 

4.55 

1.10 

3.19, 

0.67 

1.40 

11.32 

1.10, 



+0.09 
0.29 

+0.05 
0.62 

+0.05 
0.75 

+0.14 
0.62 

+0.16 
0.50 
0.03 

-0.68 

+ 1.07 
1.56 

+0.24 

-0.72 

+0.08 
0.40 

+0.03 
0.50 
1.16 
1.49 

+0.03 
0.43 

+0.13 

0.92 

+0.03 
0.32 
0.07 



(3.09: 



-1.99 >, 
-0.38^ 



(1-911?,:??) 



0.60 
+2.45 
2.35 

+0.07 

'0.33 



(9.2it?:^^) 



''PL is a power-law and NSA is a single temperature NS hydrogen atmosphere model. For 47 Tuc J and O, the PL index 
has been fixed at the nominal value obtained for the BB+PL fit as the fit parameters cannot be reliably constrained for the 
NSA+PL with the PL index as a free parameter. 

'^-Rcff calculated assuming a distance of 4.5 kpc. 

'^Spectrum may be contaminated by background/neighboring source(s). 



component fit represent averages, weighted based on the 
relative temperatures and areas of the two emitting re- 
gions. 
For the fainter MSPs, such as 47 Tuc C, M, N, Q, 



and T, a pure power-law model, with a fairly steep spec- 
trum (r > 2.5), also yields an acceptable fit. However, 
we believe that this likely stems from the fact that a 
two-temperature spectrum appears power-law-like in the 
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Fig. 3. — Representative X-ray spectra and best fit single-temperature thermal models {solid line) for the 47 Tuc MSPs. The lower panel 
for each MSP shows the best fit residuals. See the electronic edition of The Astrophysical Journal for the spectra of all 47 Tuc MSPs. 



energy range of interest (see e.g. Fig. 2 in IZavlin et al.l 
1200 21 . We also note that a broken power-law model, with 
a break at i?b « 1.0 keV, is also a good description of the 
observed emission from the 47 Tuc MSPs. In principle, 
such a spectrum could arise due to a deficit of radiating 
high-energy particles in the pulsar magnctosphere. A 
break in the spectrum is observed in young normal pul- 
sars though at m uch higher energies (a few GeV, see e.g., 
iThompsonl 11996 ). However, evidence against a purely 
non-thermal model comes from optical observations of 
PSR J0030-H0451, an MSP with a qualitatively similar 
spectrum to those of the 47 Tuc pulsars, which show that 
such a model grossly overestimates (by a factor of ^^500) 
the optic al flux upper limits whe n extrapolated to lower 
energies (jKontsevich et al.ll2003D . 

For the most luminous MSPs in the sample, namely 
47 Tuc J, O, and W, a composite NSA-hPL or BB+PL 
model yields an acceptable fit, with the PL component 
contributing -70%, --50%, and -75% of the total flux, 
respectively. The spectral fits, shown in Figure 4, were 
found to be acceptable for photon indices F ~ 1.0 — 1.4. 
We note that although an F-test does not indicate that a 
composite spectrum is statistically prefered over a pure 
PL (but with a steeper photon index of F ~ 1.5—2.0), the 
fainter thermal component is very likely genuine, given 
that the inferred values for Tofi and i?eff are very similar 
to those of the other 47 Tuc MSPs. The non-thermal X- 
rays from MSP W, and probably MSP J and O as well, 
most likely originate in an intrabinary shock (see §4) al- 
though for J and O at least a portion of this emission 



may originate in the pulsar magnetosphere. For 47 Tuc 
O, which is near the crowded center of the globular clus- 
ter, the PL component may also be due to a background 
X-ray source. The spectra of all other 47 Tuc MSPs may, 
in fact, contain a faint magnetospheric emission compo- 
nent, which is undetectable in our observations due to 
the limited number of counts. 

It is important to note that, in order to model the NSA 
and BB spectra of MSPs properly, one needs to take into 
account the rotation of the NS, which induces a modula- 
tion in the projected area of the polar cap(s) at the spin 
period and, in the case of a NSA, in the observed spec- 
trum as well due to the energy dependence of the limb 
darkening effects o f a NS atmosphere (see e.g. iR.omanil 
[19871 IZavIm et all ri996: Llovd 200|). Therefore, the 
derived temperatures and radii represent time-averaged 
values integrated over the rotation period of the pulsar. 
As a consequence, the spectral fits may underestimate 
the true flux by up to a factor of ^^2, even if we consider 
the effects of gravitational r edshift and bending of light 
near t he NS surf ace (see e.g.'Misner. Thorne. & WheeL 
1970t iPechenick . Ftaclas. & Cohen 1983; Bcloborodo; 
2002J) . Unfortunately, in order to account for these ef- 
fects, one requires prior knowledge of the relative orien- 
tation of the pulsar spin and magnetic axes and the line 
of sight, which is unavailable. 

4. VARIABILITY ANALYSIS 

Although the limited time resolution of the ACIS-S 
observations (3.2 seconds) precludes a timing analysis at 
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Fig. 4. — The X-ray spectra of 47 Tuc J, O, and W, which 
exhibit non-thermal emission. The solid line shows the best-fit 
model consisting of a power-law (dashed line) and a thermal (dot- 
dashed line) component. The lower panel for each MSP shows the 
best fit residuals. 



the millisecond level, these observations in conjunction 
with the 2000 ACIS-I observations, permit us to investi- 
gate the temporal behavior of the MSP emission over a 
large range of timescales from hours to years. 

As expected, for the bulk of the sample we find no 
statistically significant X-ray variability. For MSP O, 
however, there seems to be a gradual decrease in the 
photon count rate, by a factor of ~2 in total, especially 
at medium energies (0.6-1.5 keV), from the second to the 
fourth ACIS-S observation. There is no known physical 
mechanism which could account for such a substantial 
variation of the flux on timescales of days from an old 
MSP such as 47 Tuc O. On the other hand, such behavior 
is typical of cataclysmic variables and chromospherically 
active binaries with Lx ^ io3o-3i gj.gg g-i_ Thus, we 



conclude that this variation is most likely due to blend- 
ing of MSP O with such an X-ray source. This seems 
likely as O is at the heart of the cluster (see Fig. 1), 
where the number density of X-ray sources is large. The 
MSP R is found to be marginally variable between obser- 
vations, while MSP U exhibits possible va riability within 
the third observation IjHeinke et al.l2005|) . In both cases, 
the apparent variability is likely spurious considering the 
probability of spurious detection and the number (19) of 
trials. 

For the 12 binary MSPs, we can also examine whether 
there exist any variations in the X-ray flux as a func- 
tion of orbital position. In this sense, 47 Tuc J, O, 
R, and W are of particular interest as they are in very 
compact binaries with periods of 2.9, 3.3, 1.6, and 3.2 
hours, respectively. The former three MSPs are bound 
to very low mass degenerate companions {rric ~ 0.03 
Mq) and are eclipsed regularl y at radio wavelengths 
for 1 0-25% of their orbits (C amilo et all |2000; Frciri 
I2005|) . MSP W, on the other hand, has a significantly 
more massive (nir ^, .13 M©) main-sequence secondary 
IjEdmonds et alJl2002() and undergoes e clipses for about 
35-40% of the entire orbit l|Freirdl2005J) . For this analy- 
sis, the observations for the binary MSPs were first folded 
at the binary period using the latest values for the or- 
bital periods Pi, a nd the epochs of ascending node Tq 
IjFreire et al.ll2003l Freire et al. unpublished data). The 
net counts were then grouped in phase bins sized so as to 
allow detection of any large-amplitude variations in the 
X-ray photon count rate at the eclipse phases. For the 
same purpose, we also performed Kolmogorov-Smirnov 
and Cramer- von Mises tests on the folded but unbinned 
data, using an integer number of binary orbits. 

For 47 Tuc J, O, and R we find no detectable variation 
in the X-ray count rate at any orbital phase and are only 
able to set very crude upper limits on the amplitude of 
variation of <60-80%. In contrast, 47 Tuc W exhibits 
dramatic variations in the X-ray flux (significant at the 
99.9% level) as a function of orbital phase. This behavior 
can be naturally explained by the existence of a swept- 
back shocked stream of gas formed by interaction of the 
pulsar wind with matter from the irradiated companion 
issuing through the inner Lagrange point. The proper- 
ties of this intriguing system and its possible connection 
to low mass X-ray binaries arc described in depth by 
iBoedan ov. Grindlav. & van den Berg (200^. 

5. X-RAY LUMINOSITY VERSUS PULSAR PARAMETERS 
5.1. X-ray versus Spindown Luminosity 

Rotation-powered pulsars, including MSPs, appear to 
exhibit a linear relation between their X-ray luminosity 
Lx and rotational spindown luminosity E = ATr^IPi/P^, 
where / is the NS moment of inertia typically assumed 
to be 10^^ g cm^, with Ly ~ 10~^E (Becker fc Triimper 
119971 IPossenti et al.l l200^ . However. IGrindlav et a" 
IJ2002J) have shown that the dependence of Lx on E 
for the 47 Tuc MSPs may be significantly flatter, with 
Lx (X £;0^ 

Using the present data, we can re-examine the Lx — E 
relation for the 47 Tuc MSPs. However, the measured 
values of P and, hence, E for these MSPs are signifi- 
cantly affected by acceleration in the gravitational po- 
tential of the globular cluster. In order to obtain the 
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TABLE 4 
Derived E and Lx for the 47 Tuc MSPs. 



MSP 



(xlO^* ergs s"^) 



(xlO^" ergs s-i) 



c 


o.o5l°:Ji 


D 


0-67l°iI 


E 


3-12l°o:79 


pc 


Am+_l\l 


G(+I) 


<1.72 


Hd 


- 


I 


<7.14 


J 


3.22ll:«l 


L 


1 n4+i'^'^ 


M 


- 


N 


+ 1.49 


O 


o -,9+0.85 
■^•J^^-1.13 


Q 


1 09+0.12 


R 


9 04+3.10 
2.«4_228 


3== 


9 97+0.49 


T 


l-09ta69 


U 


3.98l°|i 


W<i 


- 


Y 


4.82±!ii 



s+2.9 

-2.4 
s+2.4 
'-2.7 
3+9.8 

-7.3 
:^+6.8 

-4.2 
5+9.3 

-7.3 
40+5.4 



2.6: 

4.i 

8.9:1 
5.51 

9.3"* 



5.6 

13. 
3.8 
3.9 
9.3 
3.9 

11.0 
8.9 
2.9 
5.1 
5.9 
3.9 



+4.6 
5.0 
12.1 
9.8 

+5.2 

3.3 

+5.5 
3.4 

+8.1 
7.6 

+5.5 
3.4 
12.1 
8.6 

+9.9 
6.5 

+3.1 
2.7 

+7.8 
4.4 

+3.3 
5.9 

+5.5 

-3.4 



^Values for E derived assuming a King model for 47 Tuc 
using the following parameters: distance D=4.5 kpc, central 
velocity dispersion v^{Q) = 11.6 km s~^, core radius r^ = 
23.1", and best fit constant gas density Ue =0.1 cm~ . 

''Values calculated using Lx = A-kB^osbT'^, with R and T 
taken from Table 2, and corrected for gravitational redshift 
assuming I + Zg = 1.31. For MSP G, Lx is for G and I 
combined. 

"^Lx apportioned based on procedure described in §2. 

•^No reliable measurements of P. 



intrinsic period derivative Pi, we first determined the 
three-dimensional position of each MSP relative to the 
cluster center by considering the projected radius and 
assuming that the observed differences in dispersion mea- 
sure within the MSP sample are solely due to a spread 
in distance along the line of sight for an assumed uni- 
form intra-cluster plasma IjFreire et al.l2001b() . The clus- 
ter acceleration term for each MSP was then computed 
assuming a King model for the cluster using a central 
velocity dispersion of Wz(0) = 11.6 ki n s~^, core radius 
r£ = 23T^j_aiid distance D ^ 4.5 kpc IjMevlan fc Mavod 
Il98fil iR-eire et al.ll200lU) . Subtracting the acceleration 
term from the observed P yields the intrinsic spindown 
rate Pi, which was then used to compute E. The result- 
ing values and their uncertainties are listed in Table 4. 
In this calculation we have not included MSP H, which 
exhibits anomalous variations in P, and MSP W, which 
currently has no measured P. 

Figure 5 shows the plot of Lx versus E for the sample 
of MSPs in 47 Tuc based on the values listed in Table 4. 
The distance used to compute Lx for the 47 Tuc MSPs 
was taken to be 4.5 kpc. For the MSP pair F and S, the 
combined Lx was apportioned based on the procedure 
described in §2. In the case of G and I the individual con- 
tribution of each MSP to the observed Lx is unknown. 
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Fig. 5. — X-ray luminosity versus spin-down luminosity {E = 
47r2/Pi/p3) for the MSPs in 47 Tuc (filled symbols) and all other 
MSPs detected in X-rays (open symbols). The X-ray luminosities 
due to polar cap (thermal), magnetospheric, and shock emission 
are shown as circles, squares, and triangles, respectively. For 47 
Tuc J and O, and J0437-4715 the individual emission components 
are plotted separately. The dashed line corresponds to the best fit 
for the 47 Tuc MSPs with thermal spectra, while the dotted line 
shows the linear relation Lx ex: W~ E. The numerical labels for 
the open symbols correspond to the values listed in the first column 
of Table 5. 



Also plotted are all MSPs detected in X-rays in the field 
of the Galaxy as well as in the globular clusters M4, M28, 
M30, NGC 6397, and NGG 6752. The latest values for 
the parameters of these MSPs are summarized in Table 
5. Here, we have chosen to compute the X-ray luminosity 
Lx — AttD^Fx in a broad energy band (0.1-10 keV) in- 
stead of the ROSAT band (0.1-2.4 keV) in order to take 
into account the substantial portion of the total flux (up 
to --80%) that is present beyond 2.4 keV for the MSPs 
with hard spectra. 

In this analysis, we have classified the MSPs based 
on the three distinct types of X-ray emission: (i) non- 
thermal magnetospheric, (ii) non-thermal shock, and (iii) 
thermal polar cap emission. In the latter category we 
have also included the MSPs with undetermined spec- 
tral properties as they have very similar luminosities to 
the pulsars with thermal spectra. In addition, all ther- 
mal luminosities have been corrected for the gravitational 
redshift, assuming Zg = 0.31. For the MSPs with known 
multi-component spectra, such as PSR J0437-4715 and 
47 Tuc J and O, the emission has been decomposed into 
the individual constituents. 

The thermal luminosities of the 47 Tuc MSPs alone fol- 
low the relation log Lx = (0.24± 1.10) log £;+(24±37). 
In this fit we have taken into account the errors in 
both Lx and E listed in Table 4. As evident in Fig- 
ure 5, due to the large uncertainties, the 47 Tuc sam- 
ple appears consistent with both the linear relation and 
a much flatter trend, such as Lx oc E^-^, predicted 
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by the iHarding fc Muslimovl l)2002(l polar cap heating 
model. It is interesting to note that both field and clus- 
ter MSPs with thermal spectra occupy the same region 
of the Lx — E diagram. This suggests that in terms of 
X-ray properties the two populations are indistinguish- 
able. 

Curiously, despite the fundamentally different emission 
mechanisms, all MSPs appear to follow the linear rela- 
tion Lx oc 1Q~^E, and have an efficiency of convert- 
ing E to Lx covering a surprisingly narrow range, with 
Lx/E ~ 10""* — 10~^. This is not expected given that 
in systems such as 47 Tuc J, O, and W, B1957-I-20, and 
J1740-5340, whose emission most likely originates in an 
intrabinary shock, Lx is determined, in part, by param- 
eters unrelated to the physics of the pulsar, such as the 
binary separation, the properties of the companion star, 
and the volume of the shocked region. The small range 
in Lx/E may stem from the fact that in all cases the 
high-energy emission is, ultimately, driven by the same 
underlying particle acceleration and pair production pro- 
cesses occurring in the pulsar magnetosphere. However, 
the fact that MSPs with E < 10^^ ergs s~^ appear to 
have thermal spectra, whereas those with E ^ 10'^^ ergs 
s~^ show strong non-thermal emission, suggests a pro- 
found difference in the conditions in the pulsar magne- 
tosphere (see §6.3). 

5.2. X-ray Luminosity versus Magnetic Field Strength 

Next, we consider how Lx depends on the pul- 
sar magnetic field as it governs the physical pro- 
cesses responsible both for magnetospheric emission, 
polar cap heating, and, to some degree, shock emis- 
sion. In Figure 6a, we plot Lx versus B\c = 9.35 x 
105(P/10-3 s)-5/2(p/io-20)i/2 Q the magnetic field at 
the light cylinder (re — cP/2tt) for a simple dipole. The 
MSPs with thermal spectra are consistent with log Lx — 
(0.39 ± 0.54) log JBic + (28 ± 3). As with the Lx - E re- 
lation, the fit is not particularly constraining du e to the 
large u ncertainties. Nonetheless, as noted by Sait o et all 
(| 19971) . it is striking that B1821-24 and B1937+21, which 
exhibit very luminous non-thermal (magnetospheric) X- 
ray emission, have the highest magnetic fields at the light 
cylinder, with B\c > 10^'^ G. It may be the case that 
these MSPs do, in fact, have heated polar caps but the 
resulting emission is several orders of magnitude fainter 
than the magnetospheric emission. Alternatively, the en- 
ergetics of the systen i may not favor substantia l polar 
cap heating (see, e.g. IHarding fc Mushmovlf2002() . Pul- 
sars with B\c < 10^ G, on the other hand, appear to 
have X-ray emission that is typically either thermal or of 
indeterminate (but probably thermal) character. 

In the case of the dipole magnetic field at the sur- 
face Bsurf = 3.2 X 10i9(PP)i/2 G, the Lx - Bsurf rela- 
tion (Fig. 66) seems even less defined with a large scat- 
ter for both the complete MSP sample and the 47 Tuc 
MSPs. The thermal Lx for the 47 Tuc MSPs is best fit 
bylogLx = (0.36±1.27)Bsurf + (22±37). The apparent 
lack of correlation between Lx and Bgnrf niay either in- 
dicate that the small scale magnetic field structure near 
the surface of th e polar caps devi ates from that of a pure 
dipole (see e.g.. lZhang fc Chengii2003') or that both the 
thermal and non-thermal X-ray emission properties of 
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Fig. 6. — (a) Lx versus magnetic field strength at the Hght 
cyUnder (Bic) for the MSPs in 47 Tuc (filled symbols) and all other 
MSPs detected in X-rays (open symbols), (b) Lx versus surface 
magnetic field strength (-Bgurf ). The dashed line corresponds to the 
best fit for the 47 Tuc MSPs with thermal spectra. The various 
symbols are defined as in Figure 5. 



MSPs are determined by the conditions much higher in 
the magnetosphere (e.g., at the light cylinder). 

5.3. Uncertainties in X-ray Luminosity and Pulsar 
Parameters 

In comparing theoretical predictions and observed re- 
lations it is necessary to examine and quantify all the 
effects that act to modify the observed parameters under 
consideration. First, it is important to emphasize that 
the X-ray luminosities plotted in Figures 5 and 6 repre- 
sent observed, not true quantities. This is especially crit- 
ical for the non-thermal magnetospheric emission, which 
is subject to beaming and, hence, viewing angle effects. 
Therefore, the measured value of Lx may greatly under- 
estimate the true X-ray luminosity. As for the thermal 
emission, a major source of uncertainty is introduced by 
the lack of knowledge of the masses and radii of the sam- 
ple of MSPs. These parameters affect Lx through the 
gravitational redshift and bending of light. In particu- 
lar, the true thermal luminosity, as measured at the NS 
surface, is larger by a factor of (1 -I- Zg)^ than the lumi- 
nosity measured by a distant observer. In our analysis, 
we have corrected for this effect by assuming the canoni- 
cal values of M = 1.4 Mq and i? = 10 km. However, for 
a reasonable range of NS masses, 1.2 -2.3 Mf?^, b ased on 
the results of iLvne et all l|2004) and lNice et all (2005) . 
we find that the true Lx may differ from the value used 
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TABLE 5 

Parameters for MSPs detected in X-rays outside of 47 Tuc. 







P 


Pi 


D 


Fx (0.1-10 keV)'' 






MSP^' 


(ms) 


(10-20) 


(kpc) 


(10-13 ergs cm-2 g-i) 


References 


1 


J0030+0451 


4.86 


<1.0'= 


0.3 


3.6 


1, 2 


2 


J0218+4232 


2.32 


7.78 


2.7 


4.2 


3, 4 


3 


J0437-4715 


5.76 


1.86 


0.139'' 


12.9 (3.0) 


5, 6 


4 


J0751+1807 


3.48 


0.73 


1.1 


0.47 


7, 8 


5 


J1012+5307 


5.26 


0.97 


>0.7T^ 


1.3 


9, 8 


6 


J1024-0719 


5.16 


<0.30'= 


<0.20<' 


0.2 


10, 11, 12 


7 


J1744-1134 


4.07 


0.71 


0.357<i 


0.2 


10, 13, 12 


8 


B1937+21 


1.56 


10.6 


3.6 


3.3 


11, 14 


9 


B1957+20 


1.61 


1.15 


2.5 


0.9 


11, 15 


10 


J2124-3358 


4.93 


1.30 


0.25 


3.6 


10, 11, 12 


11 


B1620-26 {M4) 


11.07 


~0.67'= 


1.73'^ 


0.03 


16, 17, 18 


12 


J1740-5340 (NGC 6397) 


3.60 


3.90 


2.55<i 


0.24 


19, 20, 21 


13 


B1821-24 (M28) 


3.05 


162 


5.5^ 


4.6 


22, 23 


14 


J1911-6000C (NGC 6752) 


5.28 


0.22 


4.ld 


0.01 


24 


15 


J2140-2310A (M30) 


11.02 


<209'^ 


9.0'* 


<0.007 


25, 26 



References. — (1) Lommen et al. 2000; (2) Becker & Achenbach 2002; (3) Navarro et al. 1995; (4) Webb, Olive, & Barret 
2004; (5) van Straten et al. 2002; (6) Zavlin et al. 2002; (7) Lundgren, Zepka, & Cordes 1995; (8) Webb et al. 2004; (9) Lange 
et al. 2001; (10) Toscano et al. 1998; (11) Toscano et al. 1999a; (12) Becker & Triimper 1999; (13) Toscano et al. 1999b; (14) 
Cusumano et al. 2003; (15) Stappers et al. 2003; (16) Lyne et al. 1988; (17) Richer et al. 1997; (18) Bassa et al. 2004; (19) 
D'Amico et al. 2001; (20) Bassa & Stappers 2004; (21) Grindlay et al. 2002; (22) Cognard & Backer 2004; (23) Becker et al. 
2003; (24) D'Amico et al. 2002; (25) Carreta et al. 2000; (26) Ransom et al. 2004. 

^For the globular cluster MSPs the host cluster is given in parentheses. 
'^Total flux. For J0437-4715 the value in parentheses gives the non-thermal flux. 
■^Poorly determined acceleration corrections. 

''Accurate distances/limits. Other estimates are obtained primarily from dispersion measure together with the Cordes & Lazio 
(2002) electron density model and are rather uncertain (up to a factor of ~2). 



here by up to a factor of ^-^1.8. In addition, for NS radii 
in the range 9 — 16 km, the derived Lx may deviate from 
the true value by <20%. 

The parameters E, -Bsurt, and Bic also depend on 
the MSP mass and radius, as ii^ ex MR^ and B ex 
VMR'^/R^. For the same range of NS masses (1.2 - 2.3 
M0), the assumption of 1.4 Mq may result in a deviation 

of up to ~60% for E and ^25% for -Bgurf and Bic, from 
the true values. This implies that a spread in masses 
within the sample of MSPs would tend to increase the 
scatter about the derived relations in both coordinates 
(e.g. Lx and E) and perhaps result in a skewed trend. 
Conversely, for a sample of MSPs that have very similar 
masses but are systematically more or less massive than 
the canonical 1.4 Mq, the result will be an offset in the 
derived intercept, with little effect on the slope of the log- 
log relation under consideration (e.g. log Lx — log E) . As 
for i?NS, for the range 9 — 16 km we find that E, Bsurf, 
and Bic may differ by up to a factor of ~2.6. The effect 
of the NS radius on the scatter in the trends discussed 
above is uncertain due to lack of knowledge of the true 
NS equation of state (EOS). For instance, in certain theo- 
retical models of the NS EOS a small change in Mns does 
not result in an appreciable change i n i?NSi while for oth- 
ers th e opposite is the case (see e.£f.. lLattimer fc PrakashI 
l2f)?rH for details). 

For the 47 Tuc MSP, the parameters derived from P are 
subject to an additional uncertainty arising from lack of 
knowledge of the magnitude of acceleration of the MSPs 
by the gravitational potential of the globular cluster. It 



is important to note that the model-dependent values de- 
rived for Pi are sensitive to the assumptions made regard- 
ing the gravitational potential and gas density profile of 
the cluster. For simplicity, in this paper we have derived 
Pi for each MSP assuming a simple constant gas density 
profile. However, it is likely that in reality the gas dis- 
tribution of 47 Tuc is more complex. For instance, the 
nearly identical dispersion measures of the MSP s near the 
cluster center, namely 47 Tuc F, L, O, and S (Freire et alJ 
200 3). may be interpreted as being due to the presence 
of a cavity in the plasma around the cluster core carved 
out by the winds of these pulsars. Hence, obtaining ac- 
curate values of Pi may require implementation of more 
elaborate models of the gas distribution, which will be 
explored in detail elsewhere. 

Another potential source of systematic uncertainty is 
the presence of an observational selection effect. It is 
probable that so far pulsar observations have sampled 
only the most luminous and nearby MSPs. The existence 
of X-ray faint MSPs, both in globular clusters and the 
field of the Galaxy, would mean that the observed narrow 
range of Lx for the MSPs with thermal spectra is simply 
an artifact of the limited sample of very low luminosity 
MSPs. 

Given the large uncertainties discussed above, any con- 
clusions drawn for the apparent relations between Lx 
and quantities derived from the pulsar spin parameters 
should be taken with much caution, as it is not possible, 
at present, to properly account for all effects that act to 
modify the observed trends. 
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6. DISCUSSION 
6.1. X-ray Properties of MS Ps 

The spectral fits described in §3 show that the ma- 
jority of 47 Tuc MSPs exhibit emission consistent with a 
thermal model, with no indication of any magnetospheric 
emission. It is probable that the MSPs with seemingly 
purely thermal emission do, in fact, have a fainter magne- 
tospheric X-ray component (<10^° ergs s~^). However, it 
remains unclear what determines the relative strengths 
of these two X-ray production mechanisms. Although 
viewing angle dependences and beaming of the magneto- 
spheric radiation can explain why thermal emission dom- 
inates in many MSPs, it is more likely that the relative 
strength of the two components is determined by the 
poorly understood details of the magnetospheric emis- 
si on and polar ca p heating mechanisms. 

iGrindlavet 311(120021 and Cheng fc Taa^ IpOOl have 
raised the possibility that the very narrow range of ther- 
mal luminosities and small emission areas imply the ex- 
istence of a strong multipole magnetic field near the sur- 
face. However, the small range in Lx is consistent with 
a simple dipole field if one considers general relativis- 
tic effects near the surface. In particular, gravitational 
bending of light causes >75% of the NS surface to be 
visible at any given time (see e.g. iPechenick et alJll983l: 
iPavlov fc Zavhnlll997t iBeloborodovl 120021 As a conse- 
quence, the degree of modulation in the flux which arises 
due to rotation of the NS is greatly reduced. This, in 
turn, diminishes the scatter in Lx for the sample of 
MSPs that would otherwise result from different orienta- 
tions of the magnetic and spin axes relative to the line of 
sight. Thus, although entirely plausible, the existence of 
a multipole magnetic field need not be invoked to account 
for the small range in Lx- 

We note that the inferred radii of the thermally emit- 
ting areas obtained in §3 are indeed somewhat smaller 
than the expected size of a magnetic polar cap for a 
simple dipole field, Rpc = (27ri?Ns/cP)^^^i?NS (e.g., 
Rpc = 2.6 km for a 3 ms pulsar). This likely arises 
due to the fact that we have fitted most MSP spectra 
with a single temperature model, where in reality a two- 
temperature model is more appropriate, as suggested by 
the results of §3. In adittion, due to variation of the 
projected area of the polar caps, induced by the rotation 
of the star, fits to the time-integrated spectrum result in 
an underestimate of the effective emission area. However, 
even with a two temperature spectrum and a correction 
factor for the effective area, in the case of a BB model the 
total emitting region is still substantially smaller than 
the expected polar cap area. If the NS surface does, in 
fact, radiate like a blackbody, this discrepancy may be 
indicative of either non- uniform polar cap heating or de- 
viations from a pure dipole field. For instance, the polar 
cap model of Harding & Muslimov predicts that the heat- 
ing rate by positrons shoul d be highest near the rim o f 
the polar cap (cf. Fig. 7 of lHarding fc Mush movll2002i) . 
resulting in an annulus-shaped emitting region with an 
area significantly smaller than that of the entire polar 
cap. IChen. Ruderman. fc Zhul l|1998fl have argued that a 
small emission area can arise if the MSP is an aligned ro- 
tator, resulting in a so-called "sqvieezed" polar cap. The 
outer gap model of Zhang & Cheng (2003) also predicts 
that the heated area should be significantly smaller than 



the whole polar cap due to the presence of a small scale 
multipole field near the surface. Finally, this result may 
imply that a blackbody is simply not a valid approxima- 
tion of the NS surface. The discrepancy may, thus, be a 
consequence of the fact that simplified blackbody models 
result in higher derived temperatures and smaller emis- 
sion areas due to the fact that the X-ray spectra of light- 
element atmosph eres are harder than blackbody spectra 
(|Zavlinet al.l ll9961. 

The results of this study as well as recent Chan- 
dra and XMM-Newton o bservations of nearby field 
MSPs llZavlin et alJ l200l iBecker fc AschenbachI l20(il 
iWebb et al.l l2004alhD have revealed that there are no 
clear systematic differences between the X-ray proper- 
ties of globular cluster and field MSPs, which could arise, 
for instance, due to different evolutionary paths. For ex- 
ample, the nearby field MSPs such as J0437-4715 and 
J0030-f0451, are quite similar to the bulk of 47 Tuc 
MSPs in terms of X-ray luminosity and spectral proper- 
ties. Also, B1821~24 in the globular cluster M28, is akin 
to B1937-I-21, characterized by large X-ray and spindown 
luminosities {Lx > 10^^ ergs s~^, E > 10^^ ergs s~^), 
hard non-thermal spectra, and very narrow radio and X- 
ray pulse profiles (Becker & TriiniDcr 1999). These sim- 
ilarities imply that if, in fact, multiple binary exchanges 
and accretion episodes occur for globular cluster MSPs, 
they may not significantly alter the elementary emission 
properties of the pulsar. 

6.2. The Eclipsing Binary MSPs 

iBogdanov et al.l ()2005j) have shown that the binary 
MSP 47 Tuc W exhibits large-amplitude X-ray vari- 
ability at the binary period. This modulation can be 
interpreted as being due to geometric occultations of 
an X-ray emitting intrabinary shock by the ~0.15 Mq 
main-sequence companion. On the other hand, in the 
three other systems in which radio eclipses are observed, 
namely 47 Tuc J, O, and R, we detect no such variabil- 
ity. The persistence of X-ray emission during the ra- 
dio eclipses confirms that at radio frequencies J, O, and 
R are not eclipsed by their respective binary compan- 
ions but rather by a tenuous ionized envelope formed 
due to mass loss by the companion driven by the pulsar 
wind, as in the case of PSRs B1957-H20 and J2051-0827 
(lArzoumanian. Fruchter. fc Tavloilll994l: iStaopers et alJ 

The unique X-ray behavior of 47 Tuc W is most cer- 
tainly the result of the presence of a main-sequence sec- 
ondary star instead of the typical very-low mass (~0.03 
Mq) degenerate companion found in "black widow" -like 
systems such as J, O, and R. While the latter three sys- 
tems are expected end products of low-mass X-ray binary 
evolution, 47 Tuc W is likely the result of an exchange 
encounter in which the origin al very low-m ass companion 
was ejected from the binary l)Freirell2005l and references 
therein). However, even if the current companion of MSP 
W is Roche-lobe filling, the existence of an intrabinary 
shock implies that another accretion phase cannot com- 
mence as the energetic pulsar wind is constraining the 
outflow of mass from the companion (see lBogdanov et alJ 
I2OO5I for details). 

6.3. X-ray versus Radio MSP Properties 
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In §5 we have found that the fog Lx — fog E relation for 
the thermally emitting MSPs has a slope of /3 ~ 0.2±1.1. 
Although such a trend is consistent with the prediction 
of the polar cap heating model of jHardinE & Muslimov 
l)2002r . Lx oc E^-^, we cannot exclude the possibility 
that the 47 Tuc MSPs follow the hnear Lx - E relation. 

Despite the large uncertainties, for the present sample 
of MSPs detected in X-rays we still observe a general 
trend of decrease of Lx as a function of spin down lu- 
minosity and magnetic field strength. If these parame- 
ters are, in turn, determined by the age of the MSP, it is 
tempting to speculate that Figures 5 and 6 represent evo- 
lutionary diagrams. Here we consider a possible evolu- 
tionary path of a MSP described within the framework of 
the iHardine. & Muslimov (2002) electron-positron (e^) 
pair production model. For relatively young and ener- 
getfo MSPs, such as PSR B1821-24 (r ~ 3 x 10^ yr), e± 
pairs are produced via curvature radiation (CR). These 
pairs, in turn, produce the observable X-ray/7-ray and 
radio non-thermal magnetospheric emission. Although 
some polar cap heating by a backflow of energetic par- 
ticles may occur in these objects, the non-thermal com- 
ponent (^10'^^ ergs s^^) is several orders of magnitude 
brighter, making it difficult to observe any emission from 
the polar caps. As the MSP ages, it undergoes rota- 
tional spindown and perhaps magnetic field decay, which 
result in a decline in the CR pair production rate. At 
a certain stage, the magnetic field becomes too low to 
permit CR pair production so inverse Compton scatter- 
ing (ICS) of soft X-rays by accelerated electrons becomes 
an important channel for e pair production. ICS pair 
production results in substantial polar cap heating by a 
backflow of positrons from the magnetosphere. Since the 
polar cap heating rate has a weaker dependence on the 
spindown luminosity {Lh ex E^^^) than magnetospheric 
emission {Lm oc E) , the thermal emission will eventually 
become dominant over the magnetospheric emission, as 
is observed in the majority of 47 Tuc MSPs and many 
field MSPs. At much later times, the MSP is expected 
to cross the ICS pair production dea th line (see Fig. 1 
in lHarding. Muslimov. fc: Zhangl2002|) and will no longer 
be observable as a pulsar, both in the radio and X-ray. 
In the end, as the polar caps cool, the MSP will fade 
away from view at X-ray energies as well. 

7. CONCLUSION 

We have presented a spectral analysis of the X-ray 
emission for the complete sample of known MSPs in 47 
Tuc. The X-ray spectra of the majority of MSPs are well 
described by a thermal model, a BB or a NS hydrogen 
atmosphere, with T^s -- (1-3) x 10'^ K, R^e ~ 0.1-3 km, 
and Lx -^ lO^o-^i ergs s"^ (0.3-8 keV). The small ob- 
served scatter in the thermal Lx is consistent with a sim- 
ple dipole magnetic field configuration. A two temper- 
ature thermal model, reminiscent of that observed from 



some nearby field MSPs, is also consistent with the emis- 
sion from the bulk of 47 Tuc MSPs and, in some cases, 
results in an improved fit. Such a spectrum could arise 
due to non-uniform polar cap heating or lateral diffusion 
of heat accross the NS surface away from the regions di- 
rectly heated by particles from the magnetosphere. The 
MSPs J, O, and W show a prominent non-thermal com- 
ponent in their spectra. At least part of the non-thermal 
X-rays are very likely due to interaction of the relativis- 
tic particle wind of the pulsar with matter from the com- 
panion. All other 47 Tuc MSPs exhibit no evidence of a 
significant non-thermal X-ray component. 

Our re-examination of the Lx — E relation has revealed 
that regardless of the emission mechanism responsible 
for production of the observed X-rays, the conversion 

efficiency of _E to Lx appears to be, in all cases, ~10~'* — 
10-3. 

All of the 47 Tuc MSPs, with the exception of 47 Tuc 
W, display no large amplitude variability in their X-ray 
emission on timescales of days, weeks, and years. The 
MSP W exhibits dramatic variations at the binary pe- 
riod, which can be attributed to the presence of an intra- 
binary shock that is eclipsed by the secondary star. The 
lack of long term variations in the flux of most MSPs 
could potentially serve as a crit erion for c lassification of 
unidentified soft X-ray sources l)Heinke et al...2005^1 . 

Given that most of the '--^lOO MSPs presently known 
have spin parameters similar to those of the 47 Tuc 
MSPs, we expect many of them to have predominantly 
thermal X-ray spectra. Moreover, for these MSPs, gravi- 
tational bending of light near the NS surface ensures that 
at least one of the antipodal heated polar caps is always 
observable, while for a wide range of angles between the 
magnetic and spin axes and the line of sight both polar 
caps are visible simultaneously for a large portion of the 
MSP spin period (see e.g. .NoUert et al...l989^1 . This ef- 
fect invariably results in a 100% duty cycle, even if the 
magnetospheric emission beams never sweep across the 
Earth or if the radio emission is very faint or heavily 
absorbed/scattered. Therefore, many MSPs with heated 
polar caps may, in principle, be detectable at X-ray en- 
ergies even if they cannot be observed as radio pulsars. 

Finally, based on our results, we conclude that in terms 
of the characteristics of their X-ray emission, there is no 
clear separation between MSPs in the field of the Galaxy 
and in globular clusters. 
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